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Abstract 

Based on spin-charge coupled drift-diffusion equations, which are derived from kinetic equations 
for the spin-density matrix in a rigorous manner, the electric- field- induced nonequilibrium spin 
polarization is treated for a two-dimensional electron gas with both Rashba and Dresselhaus spin- 
orbit coupling. Most emphasis is put on the consideration of the field-mediated spin dynamics for 
a model with equal Rashba and Dresselhaus coupling constants, in which the spin relaxation is 
strongly suppressed. Weakly damped electric-field-induced spin excitations are identified, which 
remind of space-charge waves in crystals. 

PACS nunibers: 72.25.Dc, 72.25.Rb, 72.10.-d 
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I. INTRODUCTION 



Spin-dependent transport phenomena are of great interest to both basic research and 
device apphcations. Especially, semiconductor-based spin electronics has been the subject 
of numerous investigations. In this field, spin rather than charge is exploited for signal 
processing. The spin-orbit interaction (SOI) opens the possibility of manipulating the spin 
of carriers using purely electrical means. Unfortunately, the very same SOI has the unde- 
sired effect of causing spin relaxation due to precession in a wave- vector dependent effective 
magnetic field, which is traced back to the SOI. For a quantum well grown on a [001] 
substrate, the Dyakonov-Perel spin relaxation^ is the most dominant effect. In general, 
however, spin relaxation depends on the details of the band structure and the relevant scat- 
tering mechanisms (see Ref. [2] and references therein). For an asymmetric quantum well, 
bulk-inversion asymmetry and structure-inversion asymmetry give rise to Dresselhaus and 
Rashba spin-orbit contributions to the Hamiltonian, respectively. The interplay between 
the linear Rashba and Dresselhaus terms causes a spin relaxation anisotropy that has been 
treated in a number of theoretical works. ^"^ For a quantum well grown along the [001] di- 
rection, the main axes of the spin- relaxation-time tensor are given by [110] and [110]. Most 
interesting is the observation that under idealized conditions and when the linear Rashba 
and Dresselhaus terms have equal strength, the relaxation of spin oriented along the [110] 
axis is totally suppressed."^ In this particular case, a conserved quantity exists, which hinders 
spin randomization.^'^ This remarkable behavior of the spin-relaxation time lead to the pro- 
posal of a nonbaUistic spin-FET.^ Other studies of the combined Rashba- Dresselhaus model 
referred to the spin- and charge- Hall effect. ^"^^ 

Recently, the field received a fresh impetus by the identification of an exact SU(2) sym- 
metry of the model, when reaching the condition of equal Rashba and Dresselhaus coupling 
strength. ^^'^^ From a theoretical point of view, an equivalent system is the Dresselhaus [110] 
model. The revealed symmetry gives rise to a massless mode with infinite lifetime at nonzero 
wave vector. Qualitative features of the associated spin pattern have been experimentally 
confirmed^^ by optical techniques that probe spin relaxation rates. Furthermore, it has been 
predicted^^'^^ that coherent spatial oscillations of the spin polarization develop in such a 
system under appropriate injection conditions. 

In this paper, we extend these interesting studies by treating spin effects under the 
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influence of an applied in-plane electric field. For the combined Rashba-Dresselhaus model 
with a SOI that is linear in k, spin-charge coupled drift-diffusion equations are derived in a 
systematic manner. Special results are presented and discussed for the special model with 
equal Rashba and Dresselhaus coupling strengths. 



II. SPIN-CHARGE COUPLED DRIFT-DIFFUSION EQUATIONS 

The effective Hamiltonian of our approach 

-f^o = + a{kyO-x - kxO-y) + PiKa^ - kyGy) (1) 

includes the Rashba spin-orbit term, which is due to the inversion asymmetry of the confining 
potential of the quantum well. In addition, there is the Dresselhaus coupling, which is present 
in semiconductors lacking bulk inversion symmetry. The model Hamiltonian refers to a two- 
dimensional semiconductor nanostructurc grown along the [001] direction. In Eq. (1), fc, 
m, and cXi {i = x, y, z) denote the in-plane wave vector, the effective electron mass, and the 
usual Pauli matrices, respectively, a and /? are the strengths of the Rashba and Dresselhaus 
spin-orbit couplings. Our total Hamiltonian encompasses also contributions stemming from 
the short-range spin-independent elastic scattering on impurities and the in-plane electric 
field E = {Ex, Ey,0). Its explicit form together with related kinetic equations for the spin- 
density matrix has been published recently. It should be noted that for the combined 
Rashba-Dresselhaus model, the field-induced spin polarization depends on the orientation 
of the in-plane electric field. ^^'^^ 

The main quantity for the theoretical analysis of spin-related phenomena is the spin- 
density matrix fy{k,k',t), which is calculated from quantum-kinetic equations. Prom this 
set of equations, spin-charge coupled drift-diffusion equations are derived for the physical 
components / = Tr/ and / = Trcr/, which arc integrated over the polar angle of k (denoted 
by the bar). In the case of weak SOI {akpr/h, phpr/h <^ 1, with kp and r being the Permi 
wave vector and elastic scattering time, respectively), the following ansatz is justified^^ 

nk.q,t) = -F(<,.t) ^^^^^ . /(fc.9.i) = -F(9.t) ^^^^^ . (2) 

where n{ek) denotes the Permi function and — h?k'^/{2m). Purthermore, we introduced 
the electron density n — J dep{e)n{e), with p{e) being the density of states of the two- 
dimensional electron gas. Applying this straightforward calculational scheme^-*^, spin-charge 



coupled drift-diffusion equations are obtained, which read in spatial coordinates Vj 
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The diffusion coefficient D and the mobility satisfy the Einstein relation 
{eD /n)dn/ dsp- Other quantities, which enter Eqs. (3) and (4), are defined by 

QiTTl X TTt^ 



(3) 



F (4) 



A = 



f 1 0^ 
1 
2 



A-- 1 
— 1 SikjLu + 



2ap 



(q;2 + I32^r, 



(5) 
(6) 



Q 



\ 



( 



s 



7 



1 

1 




\ 



(7) 



p a 
-a -p 


Eijk is the totally antisymmetric tensor in three dimensions. The spin flux J^j, which gives 
the i'th component of the particle flow with spin polarization along the axis j (i, j — x, y, z), 
has the general form that is in accordance with symmetry requirements.^^ Neglecting the 
spin-charge coupling and the influence of the electric field, we obtain a result 
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in which the spin flux J^i appears as a source of nonequilibrium spin polarization. ^^'^^ The 
set of basic equations (3) and (4), which restore published results^^ in the case of vanishing 
electric field, are solved and discussed in the Fourier space with respect to spatial coordinates. 
For the sake of a better readability of the paper, these equations are summarized in the 
Appendix. 



III. RESULTS AND DISCUSSION 



A number of quite interesting electric-field effects on the spin polarization occurring 
in systems with Rashba and Dresselhaus SOI can be studied on the basis of the coupled 
spin-charge drift-diffusion Eqs. (3) and (4) [or (Al) and (A2)]. To keep the presentation 
transparent, we restrict ourselves to spatially infinite systems by omitting any boundary 
effects. 
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A. Charge current 



The current density of charge carriers is defined by the time derivative of the dipole 
moment, which can be expressed by 

d 

j(t)^-ieV^—F{K,t)\^=o- (9) 

As we are mainly interested in the charge transport along the [110] and [110] directions, the 
spatial coordinates are rotated by introducing new wave vectors k± — {k^ ± Hy) / y/2. In 
these coordinates, the charge current along the k+ direction is given by 

j+{t) = -ie-^^F{K+, = 0, t) U^=o ■ (10) 

Taking into account Eq. (Al), we immediately obtain for the Laplace transformed current 
density 

Ms) = e^iE^j + + (3)F4s) - ^(a^ - (3')^iE_F,(s), (11) 

where = (E^ + Ey) / \/2 and F± = (F^ ± Fy) / \/2. The components of the density matrix, 
which enter this equation, refer to a spatially homogeneous system. The related (fc = 0) 
drift-diffusion Eqs. (Al) and (A2) are easily solved. Inserting the solution into Eq. (11) and 
taking into account only contributions linear to the electric field, we get 

with 7 — 13/ a. As the current is aligned along the electric field direction, there is no Hall 
current, the appearance of which would require the treatment of the nonlinear field term in 
Eq. (11), which is of higher order in the spin-orbit coupling constant. In Eq. (12), the pure 
charge current e/iE^F is complemented by a spin contribution, which exhibits a Drude-like 
frequency dependence governed by the spin-relaxation time Tg that may be much larger than 
the elastic scattering time r. Consequently, the frequency dispersion of the spin-mediated 
current contribution can set in at lower frequencies as in the ordinary Drude formula. The 
spin-induced current contribution disappears, when the Rashba and Dresselhaus coupling 
constants are equal {a — P). For the Rashba model (/3 = 0), Eq. (12) resembles the results 
derived previously for small polarons.^^ 
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B. Out of plane spin polarization for a — (3 



As another application of the spin-charge drift-diffusion Eqs. (Al) and (A2), we treat the 
evolution of an initial spin lattice produced at t — along the direction 

1 1 

fzo{K'+) = -fzoS{K,+ - Ko) + -f*QS{K,+ + Ko), fzo =| fzo \ e"^ , (13) 

with Ko being a given wave vector. Restricting ourselves to the special case a — (3, we obtain 
the exact solution 

Fz{k+, s) = /^o(^+) ^2 i^2f]2 ' (14) 
with the effective Laplace variable 

a = s — i/iE^K,-^ + Dk.^ + 2/ts, (15) 

and the frequency 

= 2V2^{fiE+ + 2iDk+). (16) 

The inverse Laplace and Fourier transformations of Eq. (14) are easily calculated and we 
obtain for the asymptotic dynamics at large times 

F^(r+, t) = exp [-D{kq - 2K)H] cos [^0^+ + <^ + lJiE+{KQ - 2K)t] . (17) 

For kq — 2K, the original static spin lattice survives and retains its shape in the steady state. 
Under this excitation condition, a long lived spin pattern is expected to occur. If kq slightly 
deviates from 2K, the electric field drives a propagating spin wave, the amplitude of which 
diminishes with time exponentially. The frequency of this wave is given by iiE+{k.o — 2K). 

C. Field- induced spin accumulation and Hanle effect 

In this subsection, the electric field-induced out-of-plane spin accumulation is studied for 
a different initial spin preparation. First, we focus on the contribution, which appears in 
a homogeneous electron gas {k — 0). From Eqs. (A7) to (AlO) given in the Appendix, we 
obtain the steady-state solution 

mr ii'^{El-El) , ^ 
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which apphes for a semiconductor with different Rashba and Dresselhaus SOI constants {a ^ 
P). This field- mediated spin accumulation depends on the orientation of the electric field 
within the plane and has the character of a second-order field effect. A spin accumulation, 
which is proportional to the electric field, occurs only in the plane. The quadratic field effect 
in Eq. (18) disappears for the pure Rashba {j3 = 0) and Dresselhaus {a = 0) system as well 
as under the condition \Ex\ — \Ey\. For sufficiently high electric fields [2/ts <S {/lEy/D], 
the out of plane spin polarization approaches the constant field-independent value — 
—4:al3mrcos{2ip)F/h^, with cp being the polar angle of the electric field. 

In addition to this field contribution of a homogeneous system, we study the response to 
a permanent harmonic spin generation of the form 

G{r+,t) ^Ge'''*+'^+'-+. (19) 

Disregarding the spin-charge couphng, the drift-diffusion Eqs. (A8) to (AlO) are analytically 
solved by 

,.[K^,u;) - G ^ ^ ^^^^^ ^^^^ ^^^^ ^ , (.U) 

where the short-hand notations 

in = iu-iiJLK+E+ + DK\, T+ = 2K+{fjLE+ + 2iDK+), r_ = 2K_ijlE_, K± = (q;±/9)t^ (21) 

were used. This solution becomes more transparent for the special case a = /3, where we 
obtain 

F,{k+,u;) = + ^ zu;-zf^E+{K+±2K) + DiK^±2Kr. (22) 

The dispersion relation of eigenmodes uj = ^Ej^{kj^ — 2K) is derived from the denominator 
A^_. This mode has the character of free carrier oscillations complemented by a spin part 
that gives rise to a soft mode at = 2K. 

Further information about the solution in Eq. (20) is obtained for a /3 and Q <^ 1/ts. 
Two limits can be distinguished in this case. Under the condition Q <^ 1/ts-, which applies 
to the steady state, Eq. (20) simplifies to 



G 

2/Ts+JilEf/D-' 



- TH^-TZ^^^ (23) 
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which can be interpreted as the electric-field analogy of the Hanle effect. ^^'^^ Another result 
is obtained in the opposite case l/r^- <S when the out-of-plane spin polarization depends 
on the orientation of the electric field 

C 

F = — (24) 

' 2/Ts + HE^ + Ey)f/{2Dy 

This solution dictates the behavior of the spin polarization in the limit a ^ f3. Both results 
agree for the special field configuration = Ey = E/\/2. 



D. Spin remagnetization waves for ot — (3 

Finally, we treat the relaxation of an initial homogeneous spin moment Fq. In this case, 
we use K = in our basic equations so that spin and charge degrees of freedom decouple 
from each other. Most interesting results are expected, when the Rashba and Dresselhaus 
SOI couplings coincide (a = 13). The set of Eqs. (A8) to (AlO) is easily solved for the 
Laplace-transformed functions. The result 

m = ^^^^^ + e-^*/- 1 cos(Q^t) - ^ sin(Q^t) } , (25) 

Fy{t) = -^^^^ + e-*/-{^^^cos(Q^t) - ^sin(Q^t)}, (26) 

F,{t) ^ e-^^/--| sin(Q^t) + F,oCOs(Q^t)|, (27) 

describes spin rotations with the frequency VIe = \/2Kjj,(Ex + Ey). Such rotations of the 
magnetic moment were studied previously both for hopping of small polarons^"^'^^ and for 
extended states in a two-dimensional electron gas.^^'^^ In analogy to space-charge waves, 
these eigenmodes are called spin-remagnetization waves. Typically, the amplitude of these 
excitations exponentially decrease with increasing time. It is a peculiarity of the special 
Rashba- Dresselhaus model that the in-plane magnetic moment is conserved: F^(t) — Fy{t) = 
Fr^Q — FyQ. This observation provides a further example for the occurrence of undamped spin 
excitations, whenever the coupling constants a and /? are equal. 



IV. SUMMARY 



Both for the Rashba-Dresselhaus model with equal coupling constants and for the Dressel- 
haus [110] model it has been recently reahzed that the relaxation of spin oriented in the [110] 
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axis is totally suppressed. As the spin lifetime becomes arbitrarily long within these models, 
the existence of a persistent spin grating has been predicted. Recent experiments^^'^^'^° in- 
deed revealed evidence supporting this interesting peculiarity of slow spin relaxation rates. 
The extreme suppression of spin relaxation in the special theoretical models for semicon- 
ductor nanostructures with SOI is due to a spin symmetry that gives rise to a soft mode at 
K — 2K. This weakly damped eigenmode leads also to interesting phenomena in a biased 
spin-orbit coupled system that was studied in this paper. Based on rigorous spin-charge 
coupled drift-diffusion equations for the components of the spin-density matrix, a number 
of field-mediated spin effects were considered: 

(i) A Drude-like spin-induced contribution was identified in the longitudinal charge cur- 
rent that disappears in the special Rashba-Dresselhaus model with equal coupling strengths. 
The frequency position and the width of the resonance are determined by the spin-relaxation 
time and the coupling constants a and (5. 

(ii) A persistent spin pattern created k — 2K is not destroyed by the electric field. 
However, under the condition of shght-off resonance, the field forces the pattern to move. 

(ill) In a homogeneous electron gas with Rashba and Dresselhaus SOI, an out-of-plane 
spin polarization develops, which has a nonlinear field character. In addition, long-lived 
remagnetization waves can be excited, whose frequency is given by a; = iiE+{k+ — 2K). 
These eigenmodes can be studied in a similar manner as space-charge waves in crystals. In 
addition, the electric-field analogy of the Hanle effect changes its character at a /3. 

(iv) The in-plane spin polarization of the special Rashba-Dresselhaus model with a — 13 is 
conserved and the spins rotate due to the electric field with the frequency VLe = \/2Kn{Ex + 

Ey). 

The experimental confirmation of the field-mediated spin effects predicted here would 
stimulate further progress both in basic research and technological innovations. 
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APPENDIX A: FOURIER TRANSFORMED DRIFT-DIFFUSION EQUATIONS 



The Fourier transformed version of our basic spin-charge coupled drift-diffusion equations 
(3) and (4) has the form 



d „ 
— - ijiEn Dk 
at 



F + iuj^-F- ^{a' - X t,E] • F) = 0, (Al) 
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with the abbreviations 
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For the Rashba model (/5 = 0), these equations agree with results derived previously.' 
Let us treat these equations in another representation characterized by the components 



K± 



Kx i Hy p, Fx i Fy 



(A6) 



V2 ' V2 

For excitations of the spin polarization that propagate exclusively along the direction 
(«;_ = 0), Eqs. (Al) and (A2) are written in the form 



d 2 
— - iijiE+K,+ + Dk^ 



F-Ua + P)k+F_ + '^{a' - p')fiE_K+F, = 0, (A7) 
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iiJ,E.K+ + Dk, H 

ot ^ Ts+ 



'2,771 

F+ + —{a + (3){iJiE+ + 2iDK+)F, 



q;-/3)^(q; + /3)//£;_F = 0, 



(A8) 



d ^ ^22' 
— - tiiE+K+ + Dk, H 

ot ^ T,_ 



'2,1/1 

F^ + —{a-(3)fiE^F, 



H — —{a + (3f{a - (3){iiE+ + 2iDK+)F = 0, 



(A9) 
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(AlO) 



with E± = {E^ ± Ey)/^/2 and 



2 _ {a + (3f 



2 



{a - 



(All) 



r,(a2 + /32)- 



For the particular case a — (3, the couphng between spin and charge degrees of freedom 
disappears and the secular equation for the spin components gives the following dispersion 
relations for eigenmodes of the biased spin system 



with K — 2ma/ti?. This result implies that there appears a field- induced undamped soft 
mode at = 2K, which reflects the presence of the spin rotation symmetry. 
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